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Recent optical conductivity <r(u) experiments have revealed an anomalous spectral distribution 
in the ferromagnetic phase of the perovskite system La\- x Sr x MnO?,. Using finite temperature 
diagonalization techniques we investigate er(u>) for a model that contains only the e s -orbital degrees 
of freedom. Due to strong correlations the orbital model appears as a generalized t-J model with 
anisotopic interactions and 3-site hopping. In the orbital t-J model <j(ui) is characterized by a 
broad incoherent spectrum with increasing intensity as temperature is lowered, and a Drude peak 
with small weight, consistent with experiment. Our calculations for two-dimensional systems, which 
may have some particular relevance for the double-layer manganites, show that the scattering from 
orbital fluctuations can explain the order of magnitude of the incoherent part of cr(a>) in the low 
temperature ferromagnetic phase. Moreover orbital correlation functions are studied and it is shown 
that x 2 -y 2 orbital order is prefered in the doped planar model at low temperature. 

PACS numbers: 75.10.-b 75.40. Cx 75.40. Gb 25.40.Fq 



I. INTRODUCTION 

The physics of doped manganese-oxide ceramics (man- 
ganites) is characterized by strong correlations and the 
complex interplay of spin-, charge-, and orbital-degrees 
of freedom as well as the coupling to the lattice, e.g. 
via Jahn- Teller coupling This complexity is di- 

rectly evident from the large number of phases in a typ- 
ical phase diagram. To quantify the different mecha- 
nisms it is helpful to analyze experiments where for cer- 
tain parameters one or the other degree of freedom is 
frozen out. Important experiments in this respect are the 
very detailed investigations of the optical conductivity of 
Lai-zSr^MnOs by Okimoto et al. |3|^]. These experi- 
ments show (a) a pseudogap in cr(uj) for temperatures 
above the Curie temperature T c (paramagnetic phase) 
with er(cj) essentially linear in u>, and (b) the evolution of 
a broad incoherent distribution in the range < lo < 1.0 
eV below T c , which still grows at temperatures below 
T c /10, where the magnetization is already close to sat- 
uration. Such a temperature dependence of o~(u>) over 
a wide energy region is quite unusual as compared with 
other strongly correlated electron systems near a metal- 
insulator phase boundary ||. Interestingly there is in 
addition a narrow Drude peak with width of about 0.02 
eV and little weight superimposed to the broad incoher- 
ent spectrum. In view of the large energy scale of 1 eV 
it is plausible that the orbital degeneracy is the source of 
this incoherent motion, since the spin degrees of freedom 
are essentially frozen out. 

Further optical studies for various 3D manganites (^-0] 
as well as for the bilayer system Lai^Sri.gM^Ov M con- 
firm the presence of the large incoherent absorption in the 
ferromagnetic state. 

Although the importance of this experiment was recog- 
nized immediately, the few theoretical studies ]9|-[TT[] were 



confined to simplified models or approximations, thereby 
ignoring important aspects of the full many body prob- 
lem. 

The aim of this work is to show that in a model which 
accounts for the orbital degeneracy, yet assumes that the 
spins are fully polarized, the broad incoherent spectral 
distribution of <j{uj), its increase with decreasing temper- 
ature, as well as the order of magnitude of o~{u>) at small 
u) can be explained. Our calculation also accounts for a 
small and narrow Drude peak as observed by Okimoto 
et al. in the range u> < 0.05 eV. This is a clear indi- 
cation of coherent motion of charge carriers with small 
spectral weight, i.e. in a model where due to the orbital 
degeneracy incoherent motion is dominant. 

The colossal magnetoresistance (CMR) of manganites 
p2| and their transport properties are usually studied on 
the basis of the double exchange (DE) Hamiltonian or the 
more general ferromagnetic Kondo lattice model (KLM) 
|l3| . The essence of these models is the high-spin con- 
figuration of e fl -electron and t2 g -cove electron spins due 
to a strong ferromagnetic Kondo exchange interaction 
K ~ leV. The kinetic energy in the partially filled e g 
band is lowered when neighboring spin are aligned lead- 
ing to a low-temperature ferromagnetic phase, while the 
high-T paramagnetic phase is disordered with a high re- 
sistivity. Although this model provides an explanation 
of the phases necessary for a qualitative understanding 
of CMR, it has been stressed that the double exchange 
mechanism is not sufficient for a quantitative description 




Furthermore we would like to stress here, that the 
KLM without orbital degeneracy also fails to explain the 
broad incoherent o~(u>) spectra observed in the ferromag- 
netic phase as T — > |||| . The KLM would instead lead 
to a sharp Drude absorption because of the ferromag- 
netic alignement [15], This underlines the importance of 
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the e g orbital degeneracy, although there are alternative 
proposals invoking strong electron-phonon coupling and 
lattice polaron formation |l|, |l4| , |il| . 

We start our discussion with a generic Hamiltonian 
for the manganite systems, i.e. the ferromagnetic Kondo 
lattice model with degenerate e g -orbitals, and derive for 
the spin-polarized case an effective model which contains 
only the orbital degrees of freedom. This orbital model 
consists of a hopping term between the same and differ- 
ent orbitals a and (3 on neighbor sites and an orbital 
interaction. Renaming a — a where a —\ or J, the 
model maps on a generalized anisotropic t-J model. The 
usual t-J model known from the cuprates appears as a 
special case. Our derivation includes the 3-site hopping 
processes, which appear as a natural consequence of the 
strong coupling limit. Although such terms do not influ- 
ence the orbital order for integer filling, they are impor- 
tant for the proper description of transport properties in 
the doped systems [[l7]-|2(J. 

Because of the complexity of the orbital model we 
present here first a numerical study of the orbital cor- 
relations and of the frequency dependent conductivity. 
Although there exist studies of the interplay of orbital 
and spin order at integer filling for LaMnC>3 plp^ ], the 
effect of doping has not been considered so far. The fi- 
nite temperature diagonalization |23|p4| ] serves here as an 
unbiased tool to study the optical conductivity and the 
change of orbital correlations as function of doping and 
temperature. 

The outline of the paper is as follows. In Section II we 
discuss the system in the frame of ferromagnetic Kondo 
lattice model with degenerate e g -orbitals and derive for 
the spin-polarized case an effective model which contains 
only the orbital degrees of freedom. In Section III we 
present first results for the optical conductivity for this 
orbital model. Apart from the frequency dependent con- 
ductivity spectra cr(u>) we also present the variation of 
the optical sum rule and the Drude weight with temper- 
ature for various doping concentrations. Furthermore we 
show that 3-site hopping processes, which appear in a 
natural way in the strong coupling limit, enhance sig- 
nificantly the Drude weight and the coherent motion of 
charge carriers. This section also contains a brief discus- 
sion of related theoretical work. The evolution of orbital 
correlations with temperature is discussed in Section IV 
for the two-dimensional model at different doping con- 
centrations. In particular we allow here for a rotation 
of the local quantization axis when calculating these cor- 
relation functions. The results presented were obtained 
by a finite temperature diagonalization method for finite 
clusters with up to N=16 lattice sites. The final Section 
contains a brief summary. 



II. MODEL 

We start from the generic ferromagnetic Kondo lattice 
model for the manganites, where the Mn e s -electrons are 
coupled to core spins S\ (formed by the t^g orbitals) and 
a local repulsion U between the electrons in the two e g 
orbitals: 

H = Hband + Hint + Hxondo- (1) 

Hband describes the hopping of e ff -electrons between sites 
i and j. We take the electrons to have two-fold orbital 
degeneracy labelled by a Roman index (a,b) and two fold 
spin degeneracy labelled by a Greek index (a, a'). Explic- 
itly, 

Hband = £ E?d\ a J iaa + ]T {ttd\ a J ib<T + H.C.). 
iaa (ij) aba 

(2) 

The hopping matrix element tf^ is a real symmetric ma- 
trix whose form depends on the choice of basis in ab space 
and the direction of the i — j bond. For the present study 
we shall use \x) ~ x 2 — y 2 and \z) ~ (3z 2 — r 2 )/^/3 as 
basis for the e g orbitals. From the Slater-Koster rules 

follows §f-f7|: \t xx \ = It; \t xz \ = &t; \t zz \ = \t in x- 
and y- direction. For a cubic lattice the matrix elements 
in z-direction are \t xx \ — \t xz \ = and \t zz \ = t. Here 
t — V 2 pcr / Aedp is determined by the Mn-d O-p hybridisa- 
tion Vdpu and the corresponding level splitting Ae^p- All 
matrix elements are negative except t xz along x-direction, 
which is positive. The level splitting AEi — E? — Ef 
can be controlled by uniaxial pressure or in the case of 
layered compounds even with hydrostatic pressure p8[ . 
We shall focus here on the orbital degenerate case, i.e. 
AE\ = 0, however we will keep the orbital energy term 
for the derivation of the orbital model. 

The large local electron-electron repulsion U is respon- 
sible for the insulating behavior for integer band-filling 

Hint = £ U a n ia ^n ial + U ab £ n ia n ib 

ia i 

+ Jab d\ aa d\ ba , diaa' d\ b a ■ (3) 
iaa' 

Here U a , U ab and J ab are the intra- and inter-orbital Hub- 
bard and exchange interactions, respectively. The rele- 
vant valences of Mn are Mn 3+ (S=2) and Mn 4 + (S=3/2) 
as was pointed out already by Zener |L3]]. The Hamilto- 
nian H banc i + Hi n t was considered by Zaanen and Olcs 
who showed that in general rather complex effective 
Hamiltonians result for transition metal ions with par- 
tially filled d shell near orbital degeneracy p9[ |. Here a 
simplified approach is preferred, with the ti g electrons 
of Mn ions forming core spins of S — 3/2, and thus we 
restrict the electron interactions in Hi nt to the e g bands. 
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Even in the case with local orbital degeneracy there is 
one lower Hubbard band which is partially filled in the 
case of hole doping. With one e g -electron per site these 
systems are Mott-insulators. Although the Jahn- Teller 
splitting can lead to a gap in the absence of U, it is not 
the primary reason for the insulating behavior for integer 



filling (30 1 . 

The interaction of itinerant electrons with the (S=3/2) 
core spins is given by 



H 



Kondo 



7' diao 



(4) 



leading to parallel alignment of the d-electron spin with 
the spin S{. Since the d-electron kinetic energy is favored 
by a parallel orientation of neighboring spins the ground 
state becomes ferromagnetic p3| . 

In the low-temperature ferromagnetic phase we may 
introduce an effective Hamiltonian which contains only 
the orbital degrees of freedom assuming a fully spin- 
polarized ferromagnet. The spin degrees of freedom can 
also be eliminated by high magnetic fields. This opens 
the possibility to investigate the orbital order indepen- 
dent of the spin degrees of freedom, and to shed light on 
the nontrivial question how the orbital order is changed 
upon doping with e g -electrons or holes as well as the pos- 
sible appearance of an orbital liquid state. 

The resulting model has similarities to the t-J model, 
where the spin-indices a and a 1 are now orbital in- 
dices. To stress the difference we use for the orbital 
indices the letters a, b and a, (3. Unitary transforma- 
tion H or i, = e~ s He s J3lj and the restriction to states 
without double occupancy leads to the orbital t-J model 
which has the following structure: 



E 



E'dtJia 



E<4 

(ij) ab 



Ijb 



H.c. 



H' 

n orb' 



(5) 

with the constraint that each site can be occupied by at 



most one electron, i.e. d\ a = d\ a {\ — nia). Here and in the 
following the index a denotes the orthogonal e g orbital 
with respect to orbital a. For one electron per site this 
model describes a Mott insulator. The orbital interaction 
H' orb follows as a consequence of the elimation of doubly 
occupied sites with energy U ~ U ab — J a b- 
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n orb — 9 j+urjj 



juu' aba(3 
1 



1 ' /-„ 



<5/3,b d~l + uadlpdjbd~s+u>a 



-$p,b d- }+ua d ib dtfd- J+VL > a 



(6) 



two types of contributions: (i) 2-site terms (u = u'), i.e. 
similar to the Heisenberg interaction in the standard t- 
J model, yet more complex because of the nonvanishing 
off-diagonal t ab , and (ii) 3-site hopping terms (u 7^ u') 
between second nearest neighbors. In the half-filled case, 
i.e. one e g -electron per site, only the 2-site interaction 
is operative and may induce some orbital order. In the 
presence of hole doping, i.e. for less than one e s -electron 
per site, both the kinetic energy and the 3-site contribu- 
tions in H' orb lead to propagation of the holes and to a 
frustration of the orbital order. 

The complexity of the model can be seen if we express 
the orbital interaction H' orb in terms of pseudospin oper- 



ators T? = |(n Jo - nib), = d\ a d ib and T ; = d^d^ 
where a and b denote an orthogonal orbital basis in the 
e s -space. To be specific we shall assume here a(b) = z(x). 
Keeping only the two-site contributions, i.e. u = u', one 
obtains an anisotropic Heisenberg Hamiltonian for the 
orbital interactions: 



5t 



H 



'(2) 



<ij> 



j-adj-bb { rp-\-rp— 
,ab 2 1 ,ba 2 \ ' 1 



■ TrTP 



t?-'t b ?(T.+T j + +T.-Tr) 

(^ a i?f - i^ a ) (TfTj + + T?Tr 

(tft\f - t^ttf) (T.+T? + TrT? 



(7) 
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b is equivalent to that stud- 
. Since wc arc interested here 



The orbital interaction H. 
ied by Ishihara et al. ^ll 
in the effect of holes on the orbital order we will gener- 
ally base our study on H' orb (6) which includes the 3-site 
hopping processes as well. 

We note that in the special case tf^ = S a bt and E? — 
equations (5) and (6) are identical to the usual t-J model 
( with J = At 2 /U and including 3-site processes). For the 
orbital model we shall adopt the same convention for the 
orbital exchange coupling, i.e. J = 4t 2 /U. 

The orbital degrees of freedom in combination with 
strong correlations (i.e. no doubly occupied sites are al- 
lowed) is expected to lead to incoherent motion of holes, 
although quasiparticle formation with reduced spectral 
weight is a plausible expectation on the basis what is 
known about the usual t-J model. 

In this paper we focus on the orbital degenerate case 
with Ef = E\ — 0. A detailed study of the influence of 
a finite level splitting will be presented elsewhere. 



III. OPTICAL CONDUCTIVITY 



Here u, u' = (±a, 0) or (0, ±a) are lattice unit vectors. 



The orbital interaction H' orb 



H 



'(2) 



H'orb consists of 



The frequency dependent conductivity consists of two 
parts p2] 
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+ 2Tre 2 D c S(cu). (8) 

The J-function contribution is proportional to the charge 
stiffness D c , which vanishes in insulators. This contribu- 
tion is broadened into a usual Drude peak in the pres- 
ence of other scattering processes like impurities which 
are not contained in the present model. The finite fre- 
quency absorption (or regular part) <j(u>) is determined 
by the current-current correlation function: 

*M = Nu) I dte^(Mt) Jx ). (9) 

For the derivation of the current operator we introduce 
twisted boundary conditions via Peierls construction 

t$ ui (A) = t£ ui cx P ^£ +U A(x)d x y (10) 

Assuming A = {A x , Ay, 0) we obtain from the kinetic 
energy operator the x-componcnt of the current operator 
j x = dH(A x )/dA x : 

3^ = ~ie E *i+u J u * 4+u«4- (11) 

j+uab 

An additional contribution follows from the 3-site term 
in Eq.(6) 

juu' aba{3 

(12) 

°ba = 5 i3,b d\ +a J\/^i+^a ~ S - b dl+uJU^dj+Wa, 

where the expression O^f (j, u, u') is an abbreviation for 
the operator within the last brackets of Eq.(6). 

We want to stress here that the orbital order in the 
undoped case is determined exclusively by the (two-site) 
orbital interaction (7), while the 3-site processes in (6) 
only contribute in the doped case. Nevertheless the mo- 
tion of holes and transport in general is influenced in a 
significant way by the 3-site current operator j x . In 
fact from studies of the t- J model it is known that cr(w) 
is not only quantitatively but even qualitatively changed 
by these terms |]l8| , ^9| . Although the 3-site terms are 
usually ignored in studies of the t-J model they are an 
important part of the strong coupling model and should 
not be droped in studies of the conductivity. One of the 
aims in our study of the orbital model is to analyse the 
effect of these 3-site terms in the model. 

A. Finite Temperature Lanczos Method (FTLM) 

For the calculation of <t(w) we use a generalization of 
the exact diagonalization technique for finite temperature 



developed by Jaklic and Prelovsek |23J . In this approach 
the trace of the thermodynamic expectation value is per- 
formed by a Monte-Carlo sampling. The current-current 
correlation function in (9): 

/>oo 

C(u;)=Re dte iut (j x (t)j x ). (13) 
Jo 

can be rewritten by introducing a complete set of basis 
functions | r) for the trace, and eigenfunctions | and 

| \Pp of H with eigenvalues EJ and Ej (here r is an 
irrelevant label, which will get its meaning below): 

R M 

CM ra \ E E e_/3£l > I I i* I ^ 

r=1 i,j'=l 

m | j x | r)S(w + El - ET) (14) 

with 

B. M 

^EE^ih^i 2 - ( 15 ) 

r— 1 i— 1 

These relation are exact if one uses complete sets. The 
FTLM is based on two approximations when evaluat- 
ing Eqs.(14) for C{u>) and (15) for the partition function 
Z: (1) The trace is performed over a restricted num- 
ber R of random states | r), and (2) the required set 
of eigenfunctions of H is generated by the Lanczos al- 
gorithm starting from the initial states | $q) =| r) and 
I ^o) = Jx I r ) I V ( r I Jx I r ) ! respectively. This proce- 
dure is truncated after M steps and yields the relevant 
intermediate states for the evaluation of C(u>). 

Detailed tests [^3| have shown that the results get very 
accurate already if R, M -C N st , where N st is the dimen- 
sion of the Hubert space. Therefore one is saving a lot of 
computational effort because only matrices of dimension 
M x M have to be diagonalized, i.e. much smaller than 
the dimension N st x N st of the full Hamiltonian H . 

Since in our problem T t z ot does not commute with the 
Hamiltonian H or b (5) the calculations for the orbital 
model are restricted to smaller clusters than in the t-J 
case, where different Sf ot subspaces can be treated sepa- 
rately. 

B. Results 

Typical results for the frequency and temperature de- 
pendence of <j{uj) for a two-dimensional 10-site cluster 
with 2 holes, i.e. corresponding to a doping concentra- 
tion x = 0.2, are shown in Fig. 1 with and Fig. 2 without 
3-site hopping terms, respectively. As characteristic fea- 
ture of the broad continuum we note: (a) the increase 
of the absorption <j{lo) with decreasing temperature, and 
(b) the width of the distribution ~ 2.5i measured 
at half- maximum. This half- width u>\/2 is independent of 
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temperature. We further note that similar calculations 
for a 3D cluster yield a slightly larger value cui/ 2 ~ 3.5i. 

If we compare this latter (3D) value with wf/f ~ 0.7 
eV found from the data of Okimoto et al. (for x = 0.175 
and T — 9K), we obtain an experimental estimate for 
the parameter t: t opt ~ 0.2 eV. For comparison a rough 
theoretical estimate based on Harrison's solid state table 
would yields t = V£ pa /Ae ~ 0.4 eV. 
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FIG. 1. Frequency dependent regular part of the optical 
conductivity o(uj) for a N=10 site plane with N e =8 electrons 
and J=0.25t (U=16t) at different temperatures. The calcu- 
lation includes the 3-site terms. The inset shows the optical 
conductivity (to(lu) including the Drude-part. The spectra 
were broadened with T=0.1t. 

Although a(uj) shows the anomalous increase with de- 
creasing temperature, we do not claim here that the or- 
bital model accounts for the full temperature dependence 
in the ferromagnetic phase (T < T c ). The description of 
the complete T-dependence certainly requires the anal- 
ysis of the degenerate-orbital KLM, i.e. the additional 
inclusion of the spin degrees of freedom. Yet we believe 
that the results can be compared with the data by Oki- 
moto et al. in the low-T limit, i.e. in the regime where 
the magnetization is saturated. 

To allow comparison with experimental data we 
present the conductivity in dimcnsionlcss form and with 
the proper dimensions of a 3D conductivity, assuming 
that the 3D system is formed by noninteracting layers. 
For a cubic lattice po = ha/e 2 . We further note that 
calculations for a small 3D cluster yield a rather similar 
<j(lu) distribution. If we consider Lai-^Sr^MnOa with 
lattice constant a = 5.5 A and h/e 2 = 4.11KJ7 we ob- 



tain po = 0.23 ■ 10 3 ficm. Hence the low frequency 
limit of the regular part of the conductivity cr(w — > 0) = 
crjf c ~ 0.3 - 0.4 • 10 3 (JJcm)- 1 . This is consistent with 
the order of magnitude for the low frequency limit of the 
incoherent part of the <j(ui) data of Okimoto et al. for 
Lai.^Sr^MnOa a l nc ~ 0.4 • 10 3 (0.3 • 10 3 ) (Qcm)" 1 for 
the doping concentrations x — 0.175 (0.3), respectively, 
at T = 10K. 

Hence we conclude that besides the energy scale also 
the absolute value of the incoherent <t(uj) spectrum is 
consistent with the orbital model. We stress that the 
value (order of magnitude) of a Q nc is essentially fixed by 
the conductivity sum rule and the scale W1/2, as long as 
the conductivity is predominantly incoherent. 
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FIG. 2. Optical conductivity as in Fig. 1 but without 
3-site contributions. Due to the stronger incoherence there is 
a smaller Drude peak (inset) and a smaller DC-conductivity. 

The insets in Fig.|l| and || show <Jq{ui) (8) with the 
Drude absorption at low frequency included, where we 
use an ad hoc chosen smearing parameter V = O.lt. This 
value corresponds to the experimental width T ~ 0.02 
eV taken from the (small) Drude peak observed by Oki- 
moto et al. Qj. While our model yields the weight of 
the Drude peak D c , it does not give T, which is due to 
extrinsic processes (e.g. scattering from impurities and 
grain boundaries). The charge stiffness D c together with 
r determines the DC-conductivity. 

Adopting the experimental estimate for T we find for 
the example in Fig.[l] (inset) for the low temperature DC- 
conductivity udc ~ 6.5 • 10 3 (Slcm) -1 and for the re- 
sistivity p ~ 0.15 • 10 -3 ficm, respectively. For com- 
parison, the experimental range of resistivities is e.g. 
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0.1 — 1.0- 10~ 3 ilcm in the ferromagnetic metallic phase of 
Lai-^Sr^MnOs |33|] . As we shall see below, the Drude 
weight and therefore the DC-conductivity depend con- 
siderably on the value of the exchange coupling J, and 
whether the 3-site hopping processes in the model are 
taken into account or not. These terms have a strong 
effect on the coherent motion of charge carriers. 
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FIG. 3. Temperature dependence of the kinetic energy S 
(solid line), the incoherent spectral weight 5*^ (dotted) and 
the Drude weights D c (dash-dotted line) and Dtcohn (dia- 
mond) for a N=10 site plane with N e =S electrons and J=0.25t 
(U=16t) in units of t/poe 2 . 

The charge stiffness D c can be determined in two ways: 
(a) using Kohn's relation |34|, or (b) via the optical sum 
rule which relates the integrated spectral weight of the 
real part of the conductivity to the average kinetic en- 
ergy: 



7TP 

a Q {u)<L> = - — (iO- 



(16) 



Here H km contains two contributions: (a) the usual ki- 
netic energy ~ t (5) and (b) the 3-site hopping processes 
~ t 2 /U in Eq. (6). The case with off-diagonal hopping 
considered here is a generalization of the sum rules for 
the t-J model ||| and for the t-J model including 3-site 
terms 18 2t|. Together with (8) this implies 



D r = - 



1 1 f°° 

— (H k x T) - — / <j(to)dcu. 
2N ire 2 J 0+ 



(17) 



In the following we shall use D c = S — S u as abbrevi- 
ation for this equation, where S denotes the sum rule 



expression and S u the finite frequency absorption. Fig- 
ure 3 shows the temperature variation of the sum rule 
(kinetic energy), the finite frequency absorption S u and 
the Drude weight D c for the model including the 3-site 
terms. 

Whereas at high temperature the sum rule is essen- 
tially exhausted by the finite frequency absorption cr(u>), 
we find at low temperatures a significant increase of the 
Drude weight. In this particular case D c contributes 
about 40% to the sum rule at T/t = 0.3. 

The conductivity data shown in Figs. 1 and 2 is char- 
acteristic for a 2D orbital liquid state, which is stabilized 
by thermal fluctuations. The significant increase of co- 
herency in Fig. 3 below T/t — 0.3 is due to x 2 -y 2 orbital 
ordering in the planar model (see discussion in Section 
IV and in particular Fig. 10). This is characteristic for 
the 2D version of the orbital t-J model and does proba- 
bly not occur in cubic systems for which an orbital liquid 
ground state was proposed [jioj] . Therefore we consider 
our T/t — 0.3 data for D c to be more appropriate for a 
comparison with the low-temperature 3D data than the 
T = values, which are enhanced due to orbital order. 

We have also calculated D c for T — using Kohn's 
relation |3~1 
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Kohr. 



1 d 2 E ($) 



2N d^> 2 



(18) 



<J>=<E> 



which yields a consistent value (Fig.||). Here $ = eaA a /h 
is the Peierls phase induced by the applied vector poten- 
tial with component A a . For the evaluation of Eq.|l8| 
we follow Ref. |l8[ ] and first search for the minimum of 
the ground state energy Eq with respect to the vector 
potential prj|] , and then calculate the curvature. Since 
this procedure is quite cumbersome, we have determined 
most of our D c data via the sum rule. The resulting 
sets of data show systematic trends as function of J and 
doping. 

From these results we expect in analogy with the t-J 
model that the single-particle electron Green's function 
is characterized by a pronounced quasiparticle peak to 
explain such a large fraction of coherent transport. 

From Fig]| we a ^ so see ^ na ^ ^ ne weight under a(u), 
i.e. S u , does not further increase for temperatures below 
T = 0.5t, although the temperature variation in Fig.l 
seems to suggest a significant further increase towards 
lower temperatures. Conductivity data for T < 0.5t is 
not shown in Figs. 1 and 2 because it shows pronounced 
discrete level structure, and probably requires larger clus- 
ters for a careful study. Nevertheless it appears that <t(uj) 
develops a pseudogap in the orbital ordered phase fl37| . 
Integrated quantities on the other hand are much less 
influenced by such effects. In view of the significant vari- 
ation of the <t(uj) continuum in the experimental data Q] 
in the range < T < T c , this suggests: (a) that the 
spectral change below T c is a combined effect of orbital- 
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and spin-degrees of freedom and hence requires the full 
Kondo lattice model (1) for its description, and (b) only 
the experimental data for T — > can be compared di- 
rectly with the orbital model. 

In the following we wish to shed more light on the role 
of the 3-site processes in the orbital Hamiltonian and in 
the current operator. The importance of the 3-site term 
becomes particularly clear from the J-dependence of the 
sum rule S (Fig^) and the charge stiffness D c (Fig.||) 
taken at temperature T/t — 0.3. These low temperature 
values of S and D c are approximatively independent of 
the exchange parameter J in the absence of 3-site terms. 
A similar observation was made earlier for the t-J model 
p8|^8[ . Moreover the sum rule S is proportional to the 
doping concentration x for the 3 cases shown (a;=0.125, 
0.2, and 0.3). When 3-site terms are taken into account 
both S and D c acquire a component which increases lin- 
early with J. These general features are fully consistent 
with results for the generalized t-J model |l8| . In partic- 
ular we see that the 3-site hopping terms introduce more 
coherence, for example for J = 0.5 D c ~ 0.012 without 
and 0.08 with 3-site terms. 
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FIG. 4. Kinetic energy S versus J for a N=10 site 2D 
cluster with N e =7 and 8 electrons, with (full) and without 
(open symbols) the 3-site terms {T/t — 0.3), respectively. 
The dash-dotted line (circles) represents a 16 site cluster with 
JV e = 14 and without 3-site terms. 

The change due to the 3-site terms is particularly large 
for the charge stiffness, which defines the weight of the 
low frequency Drude peak. The relative weight in the 
Drude peak D c /S is shown in Fig. 6 as function of J — 



At 2 /U . Small changes in J lead to considerable changes 
in the Drude weight and the coherent motion of carriers. 
For J = 0.25 (0.5) the 3-site terms lead to an increase of 
D c by a factor of 3 (5), respectively, for the 8 electron 
case. 
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FIG. 5. Drude weight D c versus J calculated via the 
sum-rule (17) for a N=10 site planar cluster with N e =7 and 8 
e g electrons, with and without the 3-site terms, respectively. 



32 16 



U/t 
4 



1.0 



0.8 



♦ ♦ N e /N=8/10 (no 3-site terms) 

❖ -O N e /N=8/10 

^-^N e /N=7/10 

A- - ^A N e /N=7/1 (no3-site terms) 

y f> 



0.6 



0.4 



i A 



0.2 



0.0 



0.0 



0.5 



1.0 

J/t 



1.5 



2.0 



7 



FIG. 6. The relative weight of the Drude peak D c /S vs. J 
with and without 3-site terms for a 10-site cluster with 2 and 3 
holes (x=0.2 and 0.3), respectively, at temperature T/t = 0.3. 

When we take the value t ~ 0.2 eV, which we deter- 
mined from the comparison of <j{u>) with the experimen- 
tal data, and U ~ 3 eV plj , we obtain U/t ~ 15 and 
J/t = At/U ~ 0.25 (J = 0.05 eV). For such a value for J 
we expect a relative Drude weight D c /S ~ 0.4, which is 
larger than the experimental ratio D c /N e ft ~ 0.2 found 
by Okimoto et al. for x — 0.175, identifying here the 
effective number of carriers N e ff j| with S. The small 
value for D c /N e ft found in the experiments suggests that 
J is not larger than the value estimated, otherwise we 
would expect a much too large relative Drude weight. A 
larger value for U would also have the effect to reduce 
D c . However, we have to stress here that certainly cal- 
culations on larger clusters must be performed, before 
possible finite size effects in D c can be quantified. 

C. Comparison with other work 

The optical conductivity for the orbital non-degenerate 
Kondo latice model was studied by Furukawa using the 
dynamical mean-field approximation Q. An important 
result of this calculation is, that the weight of intraband 
excitations within the lower exchange-split band should 
be proportional to the normalized ferromagnetic magne- 
tization M/M sat . Experimentally, however, N e ff(~ S) is 
still increasing when temperature is lowered even though 
M is already saturated. This contradiction to the pre- 
diction of the simple double-exchange model implies that 
some other large-energy-scale scattering mechanism sur- 
vives at low temperature, where the spins are frozen jjj. 
The dynamical mean field theory yields only one low en- 
ergy scale, and not two, i.e. there is no Drude peak plus 
incoherent structure. 

Recently Shiba et al. jll]] analysed the noninteract- 
ing two-band model (2) (i.e. without constraint) and 
argued that the intcrband transitions within the e g or- 
bitals may explain the anomalous absorption in the fer- 
romagnetic phase of Lai-^Sr^MnOa. While the fre- 
quency range of the interband transition is found con- 
sistent with the anomalous absorption, the structure of 
<j(lo) differs. In particular the noninteracting model leads 
to cr(u>) ~ ijj for small u>. Moreover they calculated the 
ratio S^/Dc ~ 0.85, i.e. they find more weight in the 
Drude peak than in the regular part of the conductiv- 
ity o~(uj). The ratio is found to be rather insensitive to 
doping in the range 0.175 < x < 0.3. This result differs 
considerably from the experimental data of Okimoto et 
al where this ratio is about 4 for x = 0.175. 

The nature of the orbital fluctuations in the ferromag- 
netic state and their effect on the optical conductivity 
was studied recently by Ishihara et al. ||Io|] using slave 



fermions coupled to bosonic orbital excitations. They 
deal with the concentration regime xt 3> J, where the ki- 
netic energy is expected to dominate the orbital exchange 
energy. Ishihara et al. arrive at the conclusion that the 
quasi two-dimensional nature of the orbital fluctuations 
leads to an orbital liquid in (3D) cubic systems. The 
orbital disorder is treated numerically in a static approx- 
imation. The <t(ui) spectrum obtained by this approach 
for T = O.K has a quite similar shape as the experi- 
mental spectrum at low temperature (with u>i/2 ~ 3t), 
however, because of the assumed static disorder, there is 
no Drude component. A further problem with the orbital 
liquid, as mentioned by the authors, is the large entropy 
expected for the orbital disordered state which seems to 
be in conflict with specific heat measurements Q. 

A microscopic theory of the optical conductivity based 
on a slave boson parametrization which combines strong 
correlations and orbital degeneracy was recently devel- 
oped by Kilian and Khaliullin for the ferromagnetic 
state at zero temperature. Their approach yields a highly 
incoherent spectrum due to the scattering of charge car- 
riers from dynamical orbital fluctuations, and a Drude 
peak with strongly reduced weight. However a(co) is de- 
pressed at low frequency in this calculation for the orbital 
model and an additional electron-phonon mechanism is 
invoked to obtain results similar to the experimental 
spectral distribution. This theory further accounts for 
the small values of the specific heat and therefore sup- 
ports the orbital liquid scenario for the cubic systems 

& <n> 

We note that the proposed orbital liquid state |1Cj is 

a property of the cubic systems, and not a property of 
the 2D or quasi-2D versions of the orbital model. In the 
planar model the cubic symmetry is broken from the very 
beginning and holes presumably cannot restore it. This is 
different from the physics of the t-J model, where holes 
restore the spin rotational symmetry (which is sponta- 
neously broken in the antifcrromagnetic ordered phase), 
because this symmetry is respected by the t-J Hamilto- 
nian in any dimension. 

To obtain a deeper understanding of our results for 
the conductivity, we shall analyse in the next Section 
the structure of the orbital correlations in the doped and 
undoped two-dimensional orbital model. 

IV. ORBITAL CORRELATIONS AND DOPING 
DEPENDENCE 

In the absence of holes the anisotropic Heisenberg in- 
teraction (7) will lead to an orbital order below a certain 
temperature ~ J. Doping will destroy this order and 
lead either to a disordered orbital liquid, or may gen- 
erate a new kind of ordered state which optimizes the 
kinetic energy of the holes. To investigate this question 
for an anisotropic model it is in general not sufficient 
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to calculate e.g. the correlation function < T?T? > de- 
fined with respect to the original orbital basis. In the 
following we introduce a local orthogonal transformation 
of the orbitals on different sublattices by angles <fi and ip, 
respectively. On the A-sublattice: 



1*7 
\x) 



cos(<p)\z) + sin(4>)\x) 
—sin((j))\z) + cos(4>)\x) 



This amounts to new operators, e.g. 

ft = coa^T? + «n(2^)2? 



(19) 



(20) 
(21) 



and new correlation functions, e.g. {T?T?) is defined as 



i+Rl 



co S {24,)cos{2^){T?Ti +R ) 
- Sl n(2<b) S in(2i>)(T*T* +R ) 
'Cos(2cj))sin{2ij)(T 2 z T l x +R ) 
-sin{2ct>)co S {2^){T?Ti +R ) 



(22) 



The angles (j> and tp are chosen such that the near- 
est neighbor correlation function (T?T?) takes a maxi- 
mal (ferromagnetic) value. With this convention anti- 
ferromagnetic (AF) orbital order is accounted for by the 
choice of the local quantization axis, that is by the values 
of 4> a n d ip. 

For the undoped 2D system (x-y plane) these angles are 
4> = 45° and ^ = 135° corresponding to an ^(l^) + \ z )) 
and -j^(\x) — \z)) order at low temperatures (see Fig. @). 

(T*Tt +R ) takes a value only slightly smaller than 0.25 at 
low temperature (Fig. |8|), which implies that quantum 
fluctuations are small in this state. Hence the orbital cor- 
relation function (TfT?, R ) of the undoped orbital model 
is alternating (i.e. antiferromagnetic in the pseudo-spin 
language for the orbital degrees of freedom). That the 
AF-order seems to be established at a finite temperature 
for the small clusters means that the correlation length £ 
gets larger than the systemsize L. At small temperatures 
the correlation functions have about the same value inde- 
pendent of R, i.e. the data shows no spatial decay. Above 
T/t = 0.1 the orbital correlations are small and show a 
pronounced spatial decay, which can be considered as a 
signature of a thermally disordered orbital liquid state. 

Experimentally the coexistence of antiferromagnetic 
(or staggered) orbital order and ferromagnetic spin order 
has been established in the low-doping regime of LaMnC>3 
@|||. It should be recalled that LaMn0 3 has an A- 
type antiferromagnetic spin structure JH^ , where ferro- 
magnetic layers are coupled antiferromagnetically along 
the c-axis. In the orbital ordered case the MnOg oc- 
tahedra are deformed, and Mn 3d3 X 2_ r 2 and Zd^ y 2_ r 2 
orbitals with orientation along the x- and y-axis, respec- 
tively, have been considered as relevant occupied orbitals 
H,p>|. For the planar model we find here the alternate 



occupation of — \z)) and + \z)) orbitals 

which are also oriented along the x- and ?/-axis, however 
they differ with respect to their shape (and have pro- 
nounced lobes along z) jlfj. . 
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FIG. 7. Contourplot of the rotated nearest neighbor or- 
bital correlation function (TfTt+n), R=(1,0), as function of 
the angles <j> and ip for a N=10 site planar cluster with 
N e =10 electrons, J=0.25t (U=16t) and T=0.01t. White re- 
gions correspond to positive (ferromagnetic) and black ar- 
eas to negative (antiferromagnetic) correlation functions, i.e. 
{fifi +R )>0.23 (<-0.23), respectively. The 25 contour lines 
are chosen equidistant in the intervall [-1/4, +1/4]. 
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FIG. 8. Temperature dependence of (TfT* +R ) with 
<j> = 45° and tp = 135° for nearest-neighbors and tp — 45° 
for next-nearest neighbors, respectively. Results are shown 
for different distances R for an undoped N=10 site planar 
cluster. Parameters as in Fig. 7. The strong increase for 
T/t < 0.1 is due to the onset of alternating orbital order. 

Two holes are sufficient to remove the antiferromag- 
netic orbital order and establish a ferromagnetic orbital 
order (<j> = ip = 0°) as can be seen from Figs. || and [l(]. 
This corresponds to preferential occupation of orbitals 
with symmetry x 2 — y 2 . This can be considered as a kind 
of double exchange mechanism in the orbital sector. An 
interesting feature in Fig. 10 for the doped case is the fact 
that the correlations do not show significant spatial de- 
cay even at higher temperatures where the correlations 
are small. Moreover the correlations in the doped case 
appear to be more robust against thermal fluctuations 
than in the undoped case. The characteristic tempera- 
ture (determined from the half-width in Figs.|| and |lo| ) 
is T* ~ 0.2i for the 2-hole case while for the undoped 
system T* ~ 0.05i. This trend is consistent with the 
fact that in the doped case order is induced via the ki- 
netic energy and the corresponding scale t is larger than 
J = 0.25t. Although a more detailed analysis would be 
necessary to account properly for the doping dependence. 

It has been shown recently by Ishihara et al. |^8[ that 
in the layered manganite compounds hydrostatic pressure 



leads to a stabilization of the 
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orbitals as well. 



90 




FIG. 9. Contourplot of (T?T? +R ) for R = (1,0) as func- 
tion of angles cf> and ijj for a N=10 site cluster with two holes. 
Otherwise same parameters as for Fig. 7. 
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FIG. 10. Temperature dependence of {Tif[ +R } with 
(j> — 0° and ■0 = 0° for a 10-site cluster with two e g holes. 
Parameters as in Fig. 7. The strong increase below T/t — 0.3 
signals the onset of ferromagnetic orbital order in the planar 
model. 

A comparison of the doping dependence of spatial cor- 
relations with the t-J model is given in Fig.jll] In the 
undoped case the orbital model shows almost classical 
(orbital) Neel order, whereas in the t-J model correla- 
tions are strongly reduced by quantum fluctuations. The 
T = correlation function of the t-J model shows long- 
range order (LRO) consistent with a strongly reduced 
sublattice magnetization m z ~ 0.3 (in the thermody- 
namic limit). 
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of magnitude larger than typical J values in manganite 
systems. 
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FIG. 11. Orbital correlation functions (T — > 0) as func- 
tion of distance (squared) for a 10-site and a 4 x 4 cluster 
with zero and two holes (J/t — 0.25). In the undoped case 
the correlation function is antiferromagnetic, whereas in the 
doped case orbital correlations are ferromagnetic. Antifer- 
romagnetic orbital order is indicated by full (negative) and 
open (positive) symbols, respectively. Inset: The spin corre- 
lations (SiSi +R ) for a 20-site 2D Heisenberg model and the 
t-J model with two holes are shown for comparison (t-J data 
for J — 0.4 from Ref. [47]). The two-hole case shows a rapid 
decay of correlations characteristic for an antiferromagnetic 
spin liquid. 

In the orbital model the effect of the two holes (x = 
0.125 and 0.2) is strong enough to induce ferromagnetic 
(orbital) correlations, i.e. with a prefered occupation of 
x 2 -y 2 orbitals. In the t-J model instead spin-correlations 
decay rapidly for two holes on a 20-site cluster (x — 0.1), 
which is consistent with the notion of an AF spin-liquid 
state (See inset Fig.pd|). 

The preference of ferromagnetic orbital order in the 
doped case is due to the fact that t a/3 depends on the 
orbital orientation, and t aa ^> t bb . One would expect 
that for sufficiently strong orbital exchange interaction 
J = At 2 /U one reaches a point where antiferromagnetic 
orbital interactions and ferromagnetic correlations due to 
the kinetic energy are in balance. This quantum critical 
point between FM and AF orbital order turns out to be 
at quite large J values. For the two-hole case (x = 0.2) 
we have found that this crossover happens for quite large 
orbital interaction J ~ 2, i.e. at a value about an order 



FIG. 12. Orbital occupation n x 2_ y 2 = N x 2 _ y n / 'N e of the 
x 2 -y 2 orbital as function of J in doped systems for a 10-site 
cluster with 1,2, and 3 holes. 

The relative x 2 -y 2 orbital occupation is displayed in 
Fig.|l^ for various doping concentrations as function of J 
at low temperature. This clearly shows that for realis- 
tic values for J a large occupation of the x 2 -y 2 orbital 
is obtained. The predominance of this "ferromagnetic" 
orbital order distinguishes the charge propagation in the 
orbital model from the physics in the cuprates, where the 
carriers move in an antiferromagnetic spin liquid and at 
low doping in an antiferromagnetic ordered state, respec- 
tively. Nevertheless incoherent motion persists in the or- 
bital model because of the non-diagonal hopping matrix 
elements t a ^ . 

Before closing this section, we shall explain the dif- 
ferent orbital occupancy in the doped and undoped case 
in more physical terms. The overlap of the two sets of 
orthogonal e g orbitals on neighboring sites changes with 
the angles cj) and ip. (1) In the undoped state the di- 
rect hopping between the (predominantely) occupied or- 
bitals -^(\x) + \z)) and -^(\ x ) — \ z )) on neighboring sites 
is small (Figjl3|). The hopping between these occupied 
orbitals is blocked anyhow because of Pauli's principle. 
However t ab between an occupied orbital on one site and 
an unoccupied orbital on a neighbor site is extremal. This 
leads to a maximal antiferromagnetic orbital exchange in- 
teraction, and thereby to a lowering of the energy. (2) 
In the doped case instead, the x 2 -y 2 orbital occupancy is 
prefered on all sites (Fig.[l4|), which implies a large hop- 
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ping amplitude for holes in the partially occupied band, 
whereas the off-diagonal hopping matrix element t ab and 
the orbital interaction are reduced for this choice of oc- 
cupied orbitals. Hence the main energy gain in this case 
is due to the intraband correlated motion of holes in the 
x 2 -y 2 -band. Nevertheless the off-diagonal hopping plays 
an important role in the x 2 -y 2 ordered state, and leads 
to the strong incoherent features characteristic for the 
conductivity a(u>). 




FIG. 13. Alternating -^(\x)+\z)) and -^(\x)-\z)) orbital 
order in the undoped planar system. 




FIG. 14. x 2 -y 2 orbital order in the doped phase. 

Interestingly recent experiments by Akimoto et al. 
p8| proved the existence of an A-typc antiferromag- 
netic metallic ground state in a wide doping concentra- 
tion regime of the 3D system (Lai_zNd z )i_ x Sr :E Mn03 
and in particular for Lai-^Sr^Mn O3 at high doping 
x = 0.52 — 0.58. The latter compound was investigated 
by Okimoto et al. for smaller doping x where the ground 
state is the uniform ferromagnetic state. The A-type an- 
tiferromagnet consists out of ferromagnetic layers which 



are coupled antiferromagnetically. Akimoto et al. pro- 
pose that x 2 — y 2 orbitals should be occupied in this state 
and form a pseudo-2D band, as we found it here in our 
study of the ferromagnetic planar system. An experimen- 
tal study of the optical properties in this concentration 
regime would be particularly interesting. 

V. SUMMARY AND CONCLUSIONS 

We have studied the optical conductivity, charge stiff- 
ness and the optical sum rule for an effective Hamilto- 
nian that contains only the e s -orbital degrees of freedom. 
This model follows from the more general Kondo lat- 
tice model, when the spins are ferromagnetically aligned. 
The 'orbital t-J model' is expected to describe the low- 
temperature physics of manganese oxides in the ferro- 
magnetic saturated phase. 

This work was motivated by the experimental study of 
the frequency dependent conductivity of Lai-xSr^MnOa 
by Tokura's group |||. In particular these experiments 
show in the saturated low-T ferromagnetic phase a broad 
incoherent spectrum and at small frequency in addition 
a narrow Drude peak, which contains only about 20% of 
the total spectral weight. This experiment by itself shows 
that the ferromagnetic metallic state ( although fully spin 
polarized) is highly anomalous . These features cannot 
be explained by the standard double exchange model, 
but require an additional mechanism. The most natural 
mechanism, one can think of, is the twofold degeneracy 
of the e g orbitals. 

Our results for a{u) in the 2D orbital disordered regime 
show that these characteristic features follow from the 
orbital model. In particular there is (1) a broad con- 
tinuum which decreases towards high frequency with a 
half-width W1/2 ~ 2.5 and 3.5i in two and three dimen- 
sions, respectively. (2) The absorption increases with 
decreasing temperature and (3) the absolute value of the 
incoherent part of cr(u)) at low frequency and tempera- 
ture cr™ c ~ 0-4 • 10 3 (f2cm) -1 has the correct order of 
magnitude as in Tokura's experiments. (4) Despite this 
strong incoherence there is a finite Drude peak of about 
20% of the total sum rule for small J. (5) The coherent 
motion and the Drude peak become more pronounced 
for larger values of J due to the 3-site hopping processes 
in the model. For the value J/t = 0.25 estimated here 
the relative Drude weight is 40 (20) % in the model with 
(without) 3-site hopping processes, respectively. (6) The 
low temperature values for the DC-conductivity and the 
resistivity can be estimated by assuming a width T for 
the Drude peak (taken from experiment). The additional 
scattering processes contributing to T are due to impu- 
rity or grain boundary scattering, i.e. these are extrinsic 
scattering processes which are not contained in the or- 
bital model. The important information which follows 
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from the model is the weight of the Drude peak, which 
together with T determines the DC-conductivity. 

The relatively small Drude peak in comparison with 
the incoherent part of u{ui) clearly indicates that the 
carrier motion is essentially incoherent. Nevertheless we 
expect a quasiparticle peak in the single particle Green's 
function, yet with small spectral weight. Our study of the 
frequency dependent conductivity shows that the (satu- 
rated) ferromagnetic state in the manganites has uncon- 
ventional transport properties due to the scattering from 
orbital excitations in combination with the exclusion of 
double occupancy. A detailed study of the orbital dy- 
namics in the doped system is necessary for a deeper 
understanding of these issues. 

The orbital model resembles the t-J model, which de- 
scribes the correlated motion of charge carriers in the 
cuprate superconductors. Yet an important difference 
is the cubic symmetry of the pseudospin representation 
of the orbital degrees of freedom. The hopping matrix 
elements depend on the orbital orientation, and are in 
general different in the two diagonal hopping channels. 
Moreover there is also an off-diagonal hopping matrix el- 
ement, which does not exist at all in the t-J model. At 
low temperatures doping induces (ferromagnetic) x 2 -y 2 
orbital order for realistic values for the exchange inter- 
action J in the 2D orbital model. This is in striking 
contrast to t-J physics in cuprates, where the system 
changes upon doping from antiferromagnetic long-range 
order to an antiferromagnetic spin liquid. In both models 
coherent motion coexists with strong incoherent features. 

Our calculations have shown that the orbital mecha- 
nism can explain the order of magnitude of the conductiv- 
ity at low temperature. The orbital degrees of freedom 
are certainly also important for a quantitative calcula- 
tion of the collossal magnetoresistance. This is obvious 
because in the full Kondo lattice model there is a close 
interplay between orbital and spin degrees of freedom 

1011. 
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